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I. 


Introduction 


The research conducted with the support of NASA Grant No. NGR 06-992-094 
was concerned with three main areas : 

a) Exploring the possibility of using a organic monolayer formed 
from a vapor as an insulating barrier for thin film Josephson 
junctions, 

b) Experimentally investigating the effect of an organic monola 3 r er 
on the transition temperature of a thin film superconductor , and 

c) Analyzing the geometric factors which influence Josephson 
junctions and Josephson junction interferometers. 

Most of the technical details of this research has already been 
presented in the six semi-annual progress reports. This final report 
summarizes and places in perspective these many details . Each of the three 
main research area will be discussed separately. 

II . Organic Monolayers for Josephson Barriers 

Brief summary 

Pb - organic - A1 junctions were formed which displayed Josephson 
characteristics. These junctions had approximately 1 ohm junction 
resistance at room temperature but the resistance increased with time 
when the junctions were stored in room air. The Josephson characteristics 
disappeared after aging. The cause of the aging is unknown. The organic 
layer was formed from a crucible containing stearoyl peroxide with the 
substrate heated to 175°C in a vacuum of approximately 10 ^ Torr . 

From these results- it is concluded that Josephson barriers can be 
formed from an organic vapor and that the barrier is probably only a 


monolayer thick. 
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Formation of Organic Barriers 
The original model for monolayer formation is as follows: 

A thin film of superconducting metal is first deposited on a glass 
substrate in an ultra high vacuum. The ultra high vacuum reduces the 
possibility of oxidation of the metal surface before the organic reacts 
with the surface. After the superconductor is deposited, an organic vapor 
is introduced into the vacuum. The organic is a dimer molecule made up of 
two monomer molecules as shown in the figure. 

CH„ - CH_ - CH„ - ... - M - 
3 2 2 x 

M is a reactive radical or polar group. 

X 

Monomer 

CH - CH - CH„ - . . . - M = M - . . * - CH„ - CH - CH„ 

j / Z XX Z Z j 

or 

R - M = M - R 
n x x n 

where R^ is a hydrocarbon chain of length n. 

Dimer 

The dimer must be broken up into two stable monomers which then 
react with the surface. The first layer of adsorbed monomers must be 
tightly bound to the metal surface while the succeeding layers are only 
weakly bound. This would allow heat to drive off all but the tightly 
bound fir'st layer. Thus, a monolayer of the organic would be formed. The 
thickness of the layer would depend on the length of the hydrocarbon chain. 

Basically, this model proved to be correct. However, there are a number 
of experimental difficulties which must be overcome jp order to realize 
a Josephson barrier. These are 
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a) Determination of a suitable monomer 

b) Synthesis of the dimer 

c) Breaking the dimer bond 

d) Driving off excess layers 

e) Creating the organic vapor 
Each of these will be discussed separately. 

Determination of a suitable monomer 

There are two important criteria in choosing the monomer. These 
are that the monomer strongly attach to the metal film (40-100K cal/mole) 
and that the monomer be stable, i.e. not decompose. 

In order to determine, which monomer termination would bind 
tightly to superconductors such as Pb and Sn we adapted a molecular 
orbital energy program to the study of gas-solid surface reaction. This 

v ? ~ 

work resulted in two publications which describe in detail both the cal- 
culation method and the results. In brief, the calculations showed that 
terminating groups which form an oxidative bond with surface will strongly 
bond. This lead to the following monomers: 


(I) 

CH 3 - 

CH 2 - 

• • • -CH^ ■“*" 

(II) 

CH 3 - 

ch 2 - 

... -CH 2 - Ct 

(III) 

CH 3 “ 

ch 2 - 

. . . -CII 2 - 0 - 


} 

Reactive sites 

Monomers I and III are stable while monomer II decomposes quite rapidly 
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(10 sec) to 


r _ r + c0„ 

n x n 2 




This means that in order to use monomer II, the dimer must be broken up 
on the metal film surface. 

Synthesis of the dimers 

The synthesis of the dimers for monomers I and II was relatively 
simple. The synthesis for monomer III was extremely difficult as was 
carbon 14 labeled dimers for monomer II. 

A serious attempt was made to label the No. II monomers with a 
percentage of radioactive carbon 14 atom. This was tried in order to get 
a quantitative measure of the surface coverage oh the metal films. The 
technique is relatively straight forward. The Pb is vacuum deposited and 
the organic introduced. One carbon of the organic is radioactive; giving 
off a 0.14 MeV electron. After exposing the Pb in the desired manner, 
the system is opened and the Pb sample removed and placed in a counter. 

The background count is substrated from the total count and compared with 
the count expected from a monolayer. This technique could have provided 
considerable data on the effect of the variation of deposition parameters 
on the final film. Two problems arose to prevent successful completion of 
this investigation. 1) The extreme cost of the labeled compounds and 2) 
the radiation hazzard of conducting repeated tests. 

A possible solution to the first problem is to synthesize the labeled 

compounds from simpler and cheaper compounds. This proved too time consuming 

14 

and unreliable. We abandoned the C method in favor of a trial and error 
approach to forming the monolayer. 

Synthesis of the dimer for monomer III was also quite difficult primarily 
due to decomposition of intermediate solutions and the poor mixing of the 



5 


long chain molecules. The details of the synthesis procedure are given 
in the February 1972 report. While this, dimer was successfully synthesised ^ 
the amounts were small and rather impure. 

Since dimer II proved to be usable for junction formation, we chose 
to concentrate on this more manageable organic. 

Breaking the dimer band 

Since the dimer is non-reactive , it must be broken up into two monomers. 
This can be done by heat, chemical reaction, an electron beam or by UV 
radiation. The chemical reaction method requires an additional chemical in 
the system which could cause problems. The electron beam and UV radiation 
could easily break many bonds within the monomer itself , and therefore spoil 
the monolayer formation. Heat on the other hand, is simple and effective. 

To break up the dimer two methods of heating were tried: passing 

the dimer through a heated screen and heating the substrate. The latter 
proved the most effective. 

Driving off excess layers 

Heating the substrate to between 170 and 190°C was found to drive 

off the physically absorbed second layer of organic leaving a monolayer. 

Heating also broke up the adsorbed dimers. For temperatures less 

than 170°C, e.g. 150°C, a very thick wax like build up occurs. The resistance 

of the resulting junction has a resistance greater than 100 meg ohms. This 

large resistance indicates a massive build up of organic. The measured 

capacitance of such a junction is approximately 100 pf which indicates an 

o 

insulator thickness of greater than 5000 A. If the substrate temperature 
is much above 190°C, e.g. 2I0°C, then a junction resistance of approximately 
0.08 ohms results. No capacitance null is possible with these junctions. 

These data indicate a short circuit across the junction, i.e. no insulator. 
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With the substrate temperature held between 170-190°C, junctions 
of resistance between 1 to 10 ohms can be formed. This gives indirect 
evidence that monolayers of the organic, have formed. 

Organic Vapor 

To create the organic vapor, the organic powder is heated in a ' 
quartz crucible. The temperature to which the organic is heated is not 
measured directly but the current which heats the crucible provides an 
indirect measure. The organic powder is a dimer, i.e. two monomers attached 
at their reactive ends. To form the monolayer, the dimer must be broken 
down with heat into the two monomers. Therefore, the organic must be 
sufficiently hot to cause the dimer breakdown. However , if the temperature 
becomes too hot then the organic splatters out of the crucible. A very 
narrow range of heating currents has been found suitable for proper film 
deposition. Below 27 amperes only short circuit junctions are formed. 

This indicates that dimer breakdown is not occurring. Above 31 amperes, 
spattering is a serious problem. We decided on 29 amperes. 

Josephson junction procedure 

The procedure found to best produce monolayer barriers is as follows: 

1. 1" x 1" thermal glass is cleaned and placed in a substrate holder, 

2. The vacuum system is pumped to below 2 x 10 ^ Torr - this takes 
10-15 minutes. 

3. A lead strip is vacuum deposited on the glass substrate, 

4. ' The substrate with Pb strip is then heated to approximate 180°C 

for 1.5 minutes, 

5. The substrate is exposed to organic vapor (stearoyl peroxide) 


for 20 seconds, 
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The substrate is held at 180°C for an additional 1.5 minutes, 

7. A cross strip of aluminum is vapor deposited on the organic-lead, 

8. The junction is allowed to cool and self anneal for 1-5 hours 
in the vacuum. 

III. Increasing T^ with an organic layer 
Brief summary 

A number of surface enhancement investigations which employ the 
adsorption of thin layers of gases have been reported in the literature. 

As a part of our organic monolayer studies, we observed the effects of the 
organic layer on the T of our metal films. This part of the research 
program x?as not extensively pursued but some interesting results were 
observed which indicate that the T^ of A1 was increased approximately 0.4°K. 
Concept of surface enhancement of T^ 

The principle behind surface enhanced T^ is one of increasing the net 
attractive interaction betw 7 een electrons. There are several theoretical 
models to explain or predict enhancement. One such model envisions the 
adsorbing molecule as charging or discharging the surface via electron 
transfer. This seems to explain the T enhancements and depressions caused 

c 3 

by the organic molecules adsorbed on vanadium reported by Hoffman et al. 
This model does not fit the enhancement they observed after oxidation of 
the vanadium but this model does fit the oxidation and plate changing of 

4 5, to 

Al observed by Glover and Ruhl. 

7 

Thin layers of Ge on T1 and Sn have also been observed to change 
'1 . While the sign of the change fits the above charge transfer model, 
the magnitude to the required change transfer is unrealistically large. 
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However, S, Se and Cl compounds deposited on T1 and Sn yield results 

8 

consistent with this model. 

The shift in T of Al, Sn, T1 and Bi due to the adsorption of the 
noble gases, Ar and Ne, has been explained in terms of a modification of 
the phonon spectrum. 

10-12 13 

Ginzburg and more recently Bardeen have proposed an exci ton-phonon 
interaction to increase the net attraction of the electron. In this model 
a suitable dielectric is either deposited on the metal or layers of metal 
and dielectric are formed or a thin metal is deposited on the dielectric. 

In all cases, the. dielectric serves the purpose of providing the excitons. 
While there is no direct experimental evidence to support this model, 
some of the experimental results possibly can be explained by the exciton 
model. 

In ending this section, it should be mentioned that several other 
mechanisms influence the T^ of thin superconductors. These are film 
density, disorder, crystallite size and mechanical strain. All of these 
mechanisms must be considered in interpreting experimental results. 

Observed increase in T 

c 

Increased T^ of Al was indirectly observed in the tunneling character- 
istics of Pb-Organic-Al junctions. The junctions were formed by the process 
outlined in section II. The V-I characteristics of a junction are shown in 
Fig. 1. Josephson tunneling is observed at T ^ 1.3°K. This is 0.13°K above 
the transition temperature of Al. Small bumps in the V-I curves are observed 
at 1.56°K. 



These curves show that the A1 film is definitely superconducting. But 
why is not known. In the first place, the organic was deposited on the 
Pb film and is assumed to be bonded to this layer in such a manner or to 
extract electrons from the Pb. The A1 film is deposited over this organic 
and is bonded via van der Waal forces and probably some oxidation bonding 
force, although this should not be strong. It should also be noted that 
this enhancement disappeared after storing the junction in room air for 
several days. During this time the junction resistance increased from about 
1.3 to 3.0 . The cause of this increase in junction resistance is 

not known but it is assumed to be either an increase in the tunneling 
barrier thickness or the barrier height. For thin barriers, such as those 
used here, the tunneling resistance is approximately proportional to barrier 
thickness. This would imply that the barrier would double in thickness for 
a resistance change from 1.3 to 3.0. Thus it seems more reasonable to 
conclude that changes in the barrier height (which exponentially changes the 
resistance) is the cause of the junction aging. One might also conclude that 
the barrier change is due to a surface reaction at the organic-Al interface. 
This surface reaction may then be cause of the disappearance of the enhanced 

T . 

c 

IV. Analysis of Josephson Junctions and Interferometers 
Brief Summary 

A numerical technique was developed to solve the Josephson equations 
subject to the proper boundary conditions and geometrical configuration. 
Computations on a wide variety of junctions and interferometers have provided 
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considerable insight into the functioning of these devices. In addition, 
the computer technique yields quantative data for the design of high 
sensitivity magnetometers. This includes "exact" diffraction and inter- 
ference patterns, as a function of the number of junctions, different size 
junctions and loops, and mismatch of the junctions and loops. The many 
details of these calculations are given in the fiire journal papers which 
resulted from this work. 

With the above work, it is now possible to design high sensitivity 

magnetometers and determine their sensitivity to variations in fabrication 

parameters, such as junction and loop mismatch. Thin film magnetometers 

-9 

with a sensitivity of at least 10 Gauss has been shown to be possible. 
Basic Equations 

First, the equations describing a Josephson junction will be presented. 
These will be used to formulate the equations for a Josephson junction 
interferometer. Approximate and exact equations will be presented to 
give an overall view of the theoretical work to date on Josephson junction 
interferometers . 

Josephson predicted that the current density, J(x,y), through a thin 
superconducting junction would be of the form 

J(x,y) = Jj^ sin;}) (x,y) (1) 

where J^ is the maximum current density and ^(x,y) is the pair phase across 
the junction. 

The equation 

A<Kx,y) = ~~ H x n ‘ (2) 

relates the gradient of the pair phase to the local magnetic field, d is 
the barrier thickness plus twice the London penetration depth. H, the local 
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magnetic field, is composed of the app^.ljied magnetic field plus the 
magnetic field due to the current in the superconductors. 

If one neglects the magnetic field due to the current, the total current 
through the junction as a function of "the applied magnetic field has the 
following closed form solution. 


I sin(TrO / <5 ) 

I = — ~ — j- — 4>sin (o) 

J o 


(3) 


I Q is the critical current density, J^, times the area of the junction, LW. 

$ T is the applied flux, HedL, and $ is the magnetic flux quantum 
nc/2e. <Ko), the phase in the middle of the junction, is set equal to tt / 2 ? 
consistent with maximizing the current through, the junction. 

The maximum current through the junction as a function of applied magnetic 
field is more complex if one does not neglect the self-induced magnet fields. 

To find the current Maxwell's equation 

A -IT 

A x H = — J(x,y) (4) 


must be added to the set of equations describing the junction. Combining 
Eqns. 1, 2 and 4 one obtains 


2 2 

3 4 * . 3 cp 1 . . . 

— x + —j = —x sin<f>Cx,y) 

3x 3y A 


(5) 


2 

where the .Josephson penetration depth, is equal to (nc /8iTedJ^). 

A decision must now be made as to the geometry of the junction to be solved. 

Eqn. 5 ha's only been solved for the linear and asymmetrical junction. The 
reason being that for these two geometries current only flows in the y direction, 
which is the same direction as the applied magnetic field. Thus Eqn. 5 simplifies 
to the ordinary differential equation. 


— = — \ sin<j> (y) 
dy 


( 6 ) 
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Also, the total current will now be proportional to the width of the junction 
which is in the x direction. These simplifications are not valid for 
the crossed junction. 

The boundary conditions for the linear Josephson junction are found 
using Amperes Law and Eqn. 2. 


dj) 

dy 


y=L 



y=o 


4ed 

die 


H 

e 


(7) 


= 8iredl/fic 2 w (8) 

L is the junction length. 

19 

Owen and Scalapino solved graphically Eqn. 6 subject to Eqns . 7 and 8. 

Their principal result was that the total current was proportional to length 
for junction of length less than 2A . For junctions of length greater 

J 

than 6Aj the total junction current saturated and stayed constant at 
I = This implied that the Fraunhofer diffraction pattern was a good 

approximation for . junctions of length less than 2Aj. Owen and Scalapino^ 
also predicted a linear dependence of current on magnetic field for junctions 
of length L = 10A . 

The boundary conditions for an asymmetrical junction'*’"’ are 


dj^ 

dy 


d<j> 

dy 


d(j> _ 2ed 

dy T die e 
y=L 


(9) 


dj 

dy 


y=L 


dj 

dy 


y=o 


8edfrl 

dic 2 w 


( 10 ) 


Boundary Eqns. 7, 8, 9 and 10 are important in the analysis of inter- 


ferometers . 
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The first analysis of a two junction interferometer was based on 

summing the currents through two Josephson junctions, one placed at L=o 

and the other at L=L, . The junctions were considered to be eoual 

1 1 

I sin(iT<s /<3> ) 

I = — - - Jq — (sintj) (o) + sin<j) (L^)) (11) 

J' o 

Again by neglecting the self-currents, j> (L^ ) can be solved in terms of <j>(o) 
and H by solving Eqn. 2 directly. 


ML,) - *<o) + 


( 12 ) 


Here, A is the interference loop area. Eqn. 11 then simplifies to 


I 


21 ■ sin(u? T / $ ) . 

o Jo ,eA T . . 

— j-. cos (•- — II ) 

ir 4> /$ -he e 


( 13 ) 


by assuming <j> (o) adjusts such that d>(°) + -i-- - tt/2 i.e. . current is 

maximized. In Eqn. 13 the sine term is the diffraction term and the cosine 
term is the interference term. 

DeWaele and DeBruyn Ouboter 1 ’ included in their analysis the effects 
of self-induced magnetic fields in the interference loop). To accomplish 
this, the inductance, z, of the loop times the loop current is added to 
Eqn. 12 to obtain 


cKl^) = <Ko) + A H e + 


I z 
o 


(s inrf) (o) - sintj) (L^) ) 


- ( 14 ) 


Eqn. 14 is then solved simultaneously with Eqn. 11 to obtain the interference 
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pattern. All analyses to date which include the self-induced flux have 

^ , 18,16,19 

set the diffraction term equal to one. 

The major result obtained from including the self flux in the inter- 
ference loop was that the minimums in the interference pattern do not 

go to zero current. This result is in agreement with the experimental work . 

18 

Clarke and Patterson observed that one specific geometry (the 

asymmetrical current feed interferometer) produced an interference pattern 

with increased magnetic field sensitivity. They mathematically modeled 

this case by setting sinc|)(o) = 0 in Eqn. 14. Their theoretical results 

qualitatively matched their experimental work. 

16 17 

DeWaele and DeBruyn Ouboter x * ' solved the magnetic field behavior of 
multiple point contact interferometers. They neglected junction effects, 
self-currents and diffraction effects and derived 


sin* 

r 

IN 

~2eAH 

e 

Ltfc J 


— i 

0 

[”2eAH ~j 



sin 

— 

-lie 




(15) 


16,17 

In Eqn. 15 N is the number of junctions. Their experimental results 
did not closely follow Eqn. 15, for N = 2, 3, 4, 5. 

At Colorado State University we have developed an analysis technique 
for theoretically predicting the magnetic field characteristics of multiple 
Josephson junction interferometers. The Josephson penetration depth effects 
geometry-, and junction- interference loop coupling is included in the analysis 
The mathematical formulation follows the work of Owen and Scalapino^' 
who solved the one junction case. Their results have been experimentally 
verified in detail. ’ For a symmetrical interferometer composed of N 
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linear Josephson junctions, the boundary conditions are 


d(j) 

dy 


+ 

dy 


y = L 


2N-1 


4ed 


H 


'he e 


y=o 


dtj> 

dy 


dtj) 

dy 


8TTedI/-hc^w 


y = L 


2N-1 


1 y=o 


( 16 ) 


(17) 


The coupling between junctions and interference loops is obtained from 
observing continuity of magnetic field at the junction-interference loop 
boundaries. Then Eqn. 2 can be solved directly to give the coupling betvjeen 
the phase and the first derivative of phase at the junction-interference 
loop boundaries. 


^2 ^1 + d x dy ^ L l )( - L 2 L l ) 


*4 = *3 + ^ ( V (L 4 ~ V 


(18) 


, = * x _2 dj> (2H-3) T 

*2N-2 *2N-3 ' d 1 dy L 


2N-2 L 2N-3 


The L^'s refer to the edges of the various Josephson junctions (see Fig. 1) . 
In Eqn, 18 N is the number of junctions comprising the interferometer. 

The self-induced flux in the interference loops is therefore incorporated 
into the set of equations by observing continuity of magnetic field at the 
interference loop-junction boundaries, i.e. 
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d(j) 

dy 


y=L 


1 


d <j) 
dy 


y=L 


2 


d(j> 

dy 


y=L, 


d<j) 

dy 


y =L / 


( 19 ) 


d<j> 

dy 




2N-2 


d<j) 

dy 


y=L 


2N-3 


Eqns . 6, 16, 17, 18 and 19 completely describe the critical current-magnetic 
field behavior of the interferometer. The equations can be solved numerically 
to determine the interference and diffraction patterns of symmetrical current 
feed multiple junction interferometers. 

The asymmetrical current feed interferometer has a different set of 
boundary conductions than the symmetrical current feed interferometer. They 
are 


and 


dcf) 

dy 


y=o 


2ed 

~j^7” [H e - 4 ttI/cw] 


( 20 ) 


d<j> 

dy 


2ed H 
1 e 

-he 


( 21 ) 


y=L 


2N-1 


Therefore Eqns. 6, 18, 19, 20 and 21 describe the operation of the asymmetrical 


current feed interferometer. 
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Numerical Method 

The numerical solution method is based on transforming the boundary 

val ue problem to an initial value problem. The input to the numerical 

solution method contains the estimated value of current, I, the applied 

magnetic field and the detailed dimensions of the interferometer. For 

a given applied magnetic field and current, the first derivatives of the 

phase at y = 0 and y = L are computed from Eqns. (5) and (6) or (7) 

and (8) . The value of the phase at y = 0 is iterated between zero and 2 tt 

until the first derivative of the phase at y = L r , obtained from the solution 
22 

of Eqn. (4) , equals that obtained from the solution of Eqns. (5) and (6) or 

(7) and (8). If no solution is found the selected current is greater than the 

critical current and a lower value of current is then selected. This process 

continues until a solution to the equations is found. The demarcation 

between finding a solution and not finding a solution to the equations is 

the critical current. The error in phase depends on the increment size, the 

junction size and interference loop area. A computer program flow diagram 

for a three junction interferometer is shown in Fig. 2. For each junction 

and interference loop added to the interference grating a block similar to 

the dashed block in Fig. 2 must be added. 

The computed results are presented in terms of the Josephson penetration 

depth, the length of the junctions and the area of the interference loops 

23-25 

which is in contrast with previous authors who present their interference 

patterns in terms of total current aiid inductance. 
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INPUT: INTERFEROMETER DIMENSIONS, H e , ESTIMATED I 


COMPUTE d(p(0)/dY FROM 
EQ NS. 5,6 OR 7,8 


DECREASE I 


6(0)50 


INCREASE 
.. 6 ( 0 ) - 


SOLVE = (~ f S I N 6 (V) 

FROM Y = 0 TO Y = Lj 


COMPUTE 0(L 2 ) AND d(j 
FROM EQM.-9- SOLVE - 


d 0( 


s m-fy (y) 


FROM Y=.L 2 T0 Y~L 3 


13 ©( 0 ) « 2 .\ 


HAS X BEEN 
INCREASED 


7 YES 


PRINT I, H eJ 
B{Y), 6 (Y)Y 
= 0 TO Y = L 5 , 


COMPUTE' 6(1-4) AMD d©(L 4 )/dY 
FROM EQN.9, SOLVE 

^*(17)' sW«i»(y). 

FROM Y = L 4 T0 Y-L 5 
STORE d (p £ L 5 ) / tJY ’ ■ 

DOES d6(L 5 )/clY “ d6(L 5 )/ dvl 
COMPUTED FROM EONS. 

5,6. OR 7,8 


STORE 0(Y), d<j> (Y )/dY 
FROM Y = 0 TO Y“L 5 


INCREASE I 


Figure 2 

Flow diagram for computer program 


















20 


Josephson Junctions 

The result of our work on single junctions fall into two categories: 

a) The variation of junction current with junction width and 

26 

perimeter. 

b) The boundary condition, diffraction patterns, and vortex 

/ S' 

structure of asymmetrical and crossed junctions. 

The results of these investigations are discussed separately. 

The critical current through a Josephson was reported in the literature 
to be proportional to the perimeter of the junction. This was thought 
true because of the current peaks at the edges of the length. However, 
replotting the experimental data of Pritchard and Schroen and plotting our 

re- 

calculated value (fig. 3 ) , we were able to show that the critical current 

was proportional to the width of the junction and not the perimeter. 

We have examined the boundary conditions for the asymmetrical and cross 
lb 

Josephson junctions. A striking similarity appears between the diffraction 
patterns of the two geometries. 

The boundary conditions for the asymmetrical Josephson junction 
are obtained by applying Amperes law at the boundaries. These are: 

H(o) = H - AttI/cW 
e 

H(L) = II • 
e 

These imply that no current flows in the lower part of the bottom 
superconductor nor in the upper part of the top superconductor. 

The boundary conditions for the cross junction depend on 
assumptions made about the current flow in the device. One possibility, 
which we will call case I is that part of the current flows on each side 
of the superconductor leads. The positive direction of H(o,y), R(W,Y), 
II(x,o) and Il(x,L) are assumed to be in the positive y and negative x 
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■'4 * i 



Figure 3 • 

Replotting of the data of Schroen ;md Pritcliard. 
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directions. The boundary conditions are then 


W(H + H - H (x,o)) = 4 ttI/c (22) 

W(H (x,L) - H ) =0 -s ■ - - (23) 

1 0’ i2 + H (<o,y)) = bull c (24) 

H (o,y) = 0 (25) 


L H(W,y) + WH(x,L) - L H(o,y) - WH(x,o) = 4rrl/c. (26) 
is the magnetic field due to the current flowing on the underside 
of the superconductor at y = 0. 1L ^ is the magnetic field due to 
the current flowing on the top of the superconductor at x = W. The 
important point about these boundary conditions is that there are. 
magnetic fields in the y direction in the junction. This implies that 
this current density is not uniform in the x direction. In addition 
there is no magnetic field applied in the y direction to cancel the 
self-field so the maximum critical current will be less than the 
theoretical maximum. 

A different set of boundary conditions is also applicable to the 
cross geometry if one assumes a different current path in the junction. 
The current entering the junction in the positive y direction can set 
up an image current in the negative y direction in the upper super- 
conductor. The current in the upper superconductor is then assumed to 
flow around the edge of the top superconductor at y = 0 and flow in 
the positive x direction on the top of the superconductor*. The 
boundary conditions for this current flow called case II are: 

W(H(x.L) - H )) = 0 (27) 

e 

H(o,y) = Q (28) 

L(H(W,y) + H )) = 4irl/c 


(29) 
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W(H -H(x,o))=4irI /c (30) 

e 

-WH(x,o) + L H(W,y) + WH(x,L) - L H(o,y) = 4 ttI/c. (31) 

H is the magnetic field in the negative y direction on top of the 
superconductor at x = W and is equal to 4irl/cL from Eqs. 26-31. 

Equation 29 then states that H(W,y) is zero and therefore these boundary 
conditions are identical to the asymmetrical junction boundary conditions. 

With these boundary conditions the diffraction patterns can be computed. 
These are shoT-m in Fig. 4 and Fig. 5 for L = 2A and L = 6A . 

S ■ J 
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Josephson Interferometer 
27-2 c ( 

We have used the numerical method to generate and analyze multiple 
junction interferometers (Fig. o) ■ The major results are summarized here. 

The variation of the interference patterns with change in loop area for 
three junction asymmetrical current feed interferometers are shown in Fig. 7. 
These three curves exemplify the power of the numerical method to determine 
details of the interference patterns and not just give a general outline of 
the pattern features. What is also seen is the variation of pattern amplitude 
and steepness of the slope of the main peak. The large amplitude change is 
desirable since it provides a more easily detectable signal. The steeper 
the slope the greater the sensitivity of the interferometer. For these 
3 junction interferometers a large and small amplitude mode are observed. 

This can be explained in terms of the phase and flux in the two interference 
loops . 

Typical interference patterns for 2, 3, 4, and 5 junction asymmetrical 
current feed interferometers are shown in Fig. 8. Here the total loop area 
and the junction length are held constant. Note, that the amplitude and 
steepness of the major peak increases with the number of. junction. This. is 
shown separately in Figs. 9 and 10. Junction length mismatch has been found 
not to strongly effect the 2 junction patterns. The effects of junction 
mismatch on the 3,4 and 5 junction interferometers have not been determined. 

For the two junction interferometer, the variation of pattern amplitude 
as a function of junction length is shown in Fig. 11. Here the loop area is 
a parameter. The change in interference patterns with junction mismatch 
for the 2 junctions interferometer is shown in Fig. 12. 

Along with each computed point on the diffraction and interference 
pattern curves, one also computes the current and flux distribution for the 



27 



current feed interferometer 









PER CENT CHANGE IN AMPLITUDE 


/ 


100 ; 


30 


( 5 


G 


s 


X 


25 


Ojr 




TOTAL LOOP AREA *WAjd 2 


LENGTH « 0.1 a j 


4 


£ a 

NUMBER OF JUNCTIONS 

Figure 9. Percent change in interference amplitude with an increasing 
number of junctions for asymmetrical current feed inter- 


ferometers . 




31 


O 


< 

o 


8-j 

<1 


*r“ 

III 

Q 

O 

CD 











s 


s 


s. 


s 


TOTAL LOOP AREA « 40 Xj dj 
JUNCTION LENGTH® 0.1 Aj 


c.'tegag jg yt Lg W.-.j 


£L 

6 


. ’ ’ ' NUMBER OF JUNCTIONS 

Figure 10. Maximum interference pattern slope as a function of the 
number of junctions for asymmetrical current feed inter - 


S 


<? 




ferometers . 




junction length for a 2 junction asymmetrical current feed inter 



CURRENT j/wXJ,- 







l 



Effect of unequal junction lengths on the 


interference pattern of a 2 junction symmetrical 
current feed interferometer. 



Figure 12. 



34 


entire interferometer. From all of this data, a better understanding of 
the interferometer is obtained along with a quantitative basis of designing 
and optimizing the interferometer. 



35 


V - J .* i 


References 


1. J. C. Robertson and C. W. V/ilmsen, "Extended Huckel Calculations of 
Relative Chemisorption Energies, J. Vac. Sci. and Tech. 8, 53 (1971). 

2. J. C. Robertson and C. W. Wilmsen, "Preferential Orientation of CO 
Adsorption on Ni as Determined by Extended Kuckel Calculations, J. 

Vac. Sci. and Tech., 9, 901 (1972). 

3. B. M. Hoffman, F. R. Gamble, H. M. McConnell, J. Am. Chemical Society 
89, 27 (1967). 

4. R. E. Glover, III and M. D. Sherrill, "Changes in Superconducting 
Critical Temperature Produced by Electrostatic Changing," Phys. Rev. 
Letters 5, 248 (1960). 

5. R. E. Glover, III and W. Ruhl, 10th Inti. Conf . on Low Temp. Phys. , 

Moscox? (1966) . 

6. W. Ruhl, "Zur Anderung der Uber gangs temperatur Verschiedener Supraleiter 
durch Elelctronenentzug , " Zeitshrift fur Physik 196, 464 (1966). 

7. D. G. Naugle, "The Effect of Very Thin Ge Coating on the Superconducting 
Transistion of Thin Sn and T1 Films," Phys. Letters 25A, 688 (1967). 

8. P. Hilsch and D. G. Naugle, "The Influence of Dielectric Films on the 
Superconductivity of Thin Sn and T1 Filins," Zeitschrift fur Physik 
201, 1 (1967). 

9. D. G. Naugle, J. W. Balser and R. E. Allen, "Evidence for a Surface- 
Phonon Contribution to Thin-Film Superconductivity Depression of T . 
by Noble-Gas Overlayers," Phys. Rev. B 7, 3028 (1973). 

10. V. L. Ginzburg, "On Surface Superconductivity," Phys. Letters. 13, 101 
(1964) . 

11. V. L. Ginzburg, "Concerning Surface Superconductivity," J. Exptl. Theoret. 
Phys. (U.S.S.R.) 47, 2318 (1964). 

12. V. L. Ginzburg, "The Problem of High Temperature. Superconductivity," 
Contemp. Phys. 9, 355 (1968). 

13. David Allinder, James Bray and John Bardeen, "Model for an Excition 
Meehanisai of Superconductivity," Phys. Rev. B 7, 1020 (1973). 

14. C. S. Owen and D. J. Scalapi.no, "Vortex Structure and Critical Currents 
in Josephson Junctions," Phys. Rev. 164, 538 (1967). 

15. D. L. Stuehm and C. W. Wilmsen, "Diffraction Patterns and Vortex Structure 
of Asymmetrical and Cross Josephson Junctions," J. Appl. Phys. to be 
published. 

16. A. Th. A. M. DeWaele and R. DeEruyn Ouboter, "Quantum Interference 
Phenomena in Point Contacts Between Two Superconductors," Physica 41, 

255 (1969). 


36 


17. A. Th. A. M. DeWaele and R. DeBruyn Ouboter, Progress in Low Temperature 
Physics , edited by C. J. Gorter , North-Holland Publ. Co., Amsterdam. 

18. J. Clarke and J. L. Patterson, "Josephson Junction Amplifier," Appl. 

Phys. Letters 19, 451 (1971). 

19. T. A. Fulton, L. N. Dunkleberger and R. C. Dynes, "Quantum Interference 
Properties of Double Josephson Junctions," Phys. Rev. 6, 855 (1972). 

20. Walter Schroen and J. Paul Pritchard, Jr. , "Maximum Tunneling Super- 
currents through Josephson Barriers," J. Appl. Phys. 40, 2118 (1969). 

21. Klaus Schwidtal, "Type-I and Type-II Superconductivity in Wide Josephson 
Junctions," Phys. Rev. B 2, 2526 (1970). 

22. J. M. Rowell, "Magnetic Field Dependence of the Josephson Tunnel Current," 
Phys. Rev. Letters 11, 200 (1963). 

23. J. Clarke, "Electronics with Superconducting Junctions," Phys. Today 
24, 30 (1971). 

24. 0. Doyle, "Josephson Junctions Leave the Lab. Electronics 44, 38 
(1971) . 

25. J. E. Zimmerman, Paul Thiene and J. T. Harding, "Design and Operation 
of Stable rf-Biased Superconducting Point-Contact Quantum Devices, and 
a Note on the Properties of Perfectly Clean Metal Contacts," J. Appl. 

Phys. 41, 1572 (1970). 

26. D. L. Stuehm and C. W. Wilmsen, "Geometrical Dependence of the Maximum 
DC Josephson Current," J. Appl. Phys. 42, 869 (1971). 

27. D. L. Stuehm and C. W. Wilmsen, "Three Josephson Junction Interferometer," 
Appl. Phys. Letters 20, 456 (1972). 

28. D. L. Stuehm and C. W. Wilmsen, "Critical Current Periodicity of Josephson 
Junction Interferometers," J. Appl. Phys. 43, 3588 (1972). 

29. D. L. Stuehm and C. W. Wilmsen, "Josephson Junction Interference Gratings," 
J. Appl. Phys. 44, 2881 (1973). 



